Modified Newton's gravity in Finsler Space as a 
possible alternative to dark matter hypothesis 



oo 
o 
o 



o 
cr. 



(N 



(N 

o 

oo 
o 



X 



Zhe Chang and Xin Li 

Institute of High Energy Physics, Chinese Academy of Sciences 
P. 0. Box 918(4), 100049 Beijing, China 



Abstract 



A modified Newton's gravity is obtained as the weak field approximation of 
the Einstein's equation in Finsler space. It is found that a specified Finsler struc- 
ture makes the modified Newton's gravity equivalent to the modified Newtonian 
dynamics (MOND). In the framework of Finsler geometry, the flat rotation curves 
of spiral galaxies can be deduced naturally without invoking dark matter. 
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There are a great variety of observations which show that the rotational velocity 
curves of all spiral galaxies tend to some constant valuespQ. These include the Oort 
discrepancy in the disk of the Milky Way [2], the velocity dispersions of dwarf Spheroidal 
galaxies[3], and the flat rotation curves of spiral galaxiesjl]. These facts violate sharply 
the prediction of Newtonian dynamics or Newton's gravity. 

The most widely adopted way to resolve these difficulties is the dark matter hy- 
pothesis. It is assumed that all visible stars are surrounded by massive nonluminous 
matters. Though it explains the flat rotation curves of spiral galaxies, the hypothesis 
has its own weakness. No theory predicts these matters, and they behave in such ad 
hoc way. There are a lot of possible candidates of dark matter (such as axion, neutrino 
et al), but none of them satisfactory. Up to now, all of them either undetected or 
excluded by observation. 

Because of these troubles induced by dark matter, some models have been built 
for alternative of the dark matter hypothesis. Their main ideas are to assume that 
the Newtonian gravity or Newton's dynamics is invalid in galactic scale. In particular, 
two models explain well the flat rotational curves successfully without invoking dark 
matter. One is higher-order gravitational theory [S]. The gravitational potential was 
supposed j6] of the Yukawa form, 

^(r) = --^-(l + ae^). (1) 
r(l + a) ' 

Another is the famous MOND [7]. It assumed that the Newtonian dynamics does not 
hold in galactic scale. The particular form of MOND is given as 



a oJ (2) 
lim a(x) = 1, lim ri(x) = x, 

where a is at the order of 10~ 8 cm/s 2 . At beginning, as a phenomenological model, 
MOND explains well the flat rotation curves with a simple formula and a new param- 
eter. In particular, it deduce naturally a well-known global scaling relation for spiral 
galaxies, the Tully-Fisher relation[8]. By introducing several scalar, vector and ten- 
sor fields, Bekensteinjn] rewrote the MOND into a covariant formulation. He showed 
that the MOND satisfies all four classical tests on Einstein's general relativity in Solar 
system. 

These models seem appealing in theoretical interest and fit well the empirical data 
of the flat rotation curves. However, it has been pointed out that solving dark matter 
problem by means of metric theories is too difficult. The Yukawa term ([I]) is mediated 
by a spin 1 vector particle [6]. Such a contribution (TD] may be proportional to baryonic 
charge. The baryonic charge varies from one body to another, this violates the weak 
equivalence principle. MOND, as a very successful phenomenological model, still face 
problems yet. Reproducing the MOND force law[ll] requires any completely stable, 
metric-based theory of gravity to become conformally invariant in the weak field limit, 
and the prospects for a formulation with a very weak instability. 
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In this Letter, we present another possible alternative to the dark matter hypothesis. 
The major property of the flat rotational curves is the velocity of a particle in a circular 
orbit around a finite spiral galaxy becomes independent of the radius of the orbit at 
large radii. This empirical fact makes us consider the Finsler geometry. In Finsler 
geometry, the intrinsic curve is not only the function of position but also the function 
of velocity. And in the framework of Finsler geometry, the four- velocity vector is treated 
as independent variable [T2] . The explicit DISIM{,(2) invariant Finslerian line element 
and the respective Lie algebra was first proposed thirty years ago[T3"l [TH [T5] . Finsler 
geometry is a natural and fundamental generalization of Riemann geometry. 

The gravity in Finsler space has been investigated for a long time[T6 l fT7 l fT8 l fT9]. In 
this Letter, we begin with the field equation in Berwald space-a special space in Finsler 
space. The gravitational field equation in Berwald- Finsler space has been written down 
explicitly [20], 



We will study the weak field approximation of the field equation ([3]). Results obtained 
should be compared with Newton's gravity. We wish the modified Newton's gravity 
got as the weak field approximation of the field equation (J3J) in Finsler space describes 
the odd behavior of the rotation curves of spiral galaxies. 

Before dealing with the field equation, we introduce some basic notation of Finsler 
geometry [21J. Denote by T X M the tangent space at x G M, and by TM the tangent 
bundle of M. Each element of TM has the form (x,y), where x G M and y G T X M. 
The natural projection 7r : TM — > M is given by ir(x, y) = x. A Finsler structure of 
M is a function 

F : TM -> [0,oo) 

with the following properties: 

(i) Regularity: F is C°° on the entire slit tangent bundle TM\0. 

(ii) Positive homogeneity : F(x, Xy) = XF(x,y) for all A > 0. 

(iii) Strong convexity: The n x n Hessian matrix 

9 ""- dydy { 2 

is positive-definite at every point of TM\0. 

Finsler geometry has its genesis in integrals of the form 

r „ ( 1 „ dx 1 dx n \ , 

Throughout the Letter, the lowering and raising of indices are carried out by the 
fundamental tensor g^ v defined above, and its inverse g^ u . 
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In Finsler manifold, there exist a unique linear connection - the Chern connection [22] 
It is torsion freeness and metric-compatibility, 
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V 



(5) 



where 7" is the formal Christoffel symbols of the second kind with the same form 
of Riemannian connection, Ntf is defined as N^ L = ^ a y a — A l ^ x 'j^ l3 y a y 13 and A\^ u = 
^■^-^^(F 2 ) is the Cartan tensor (regarded as a measurement of deviation from 
Riemannian Manifold). The curvature of Berwald - Finsler space is given as 
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(6) 



The Ricci tensor on Finsler manifold was first introduced by Akbar-Zadeh[23j. In 
Berwald-Finsler space, it reduces to 



Ric 



R a 1 



(7) 



It is manifestly symmetric and covariant. Apparently the Ricci tensor will reduce to 
the Riemann-Ricci tensor if the Cartan tensor vanish identically. The trace of the Ricci 
tensor gives the scalar curvature S = g^ u Ric^ u . We starts from the second Bianchi 
identities on Berwald-Finsler space 



0. 



where the | means the covariant derivative. The metric-compatibility 

flV|a = and 9^\a = °> 
and contraction of (JHJ) with gives that 

=0. 

Lowering the index a and contracting with g aA , we obtain 



Av\[i 



Ric 



1 c 



_i_ J _ p> a 

1 2 a 



Ft a 



0. 



where 



B 



fivafi 



-A, uX R e \py e /F. 



(9) 



(10) 



(12) 



Thus, the counterpart of the Einstein's field equation on Berwald - Finsler space takes 
the form 



— R a 

2 a iiv 



D 



fi va 



8ttGT, 



(13) 



To get a modified Newton's gravity, we consider a particle moving slowly in a weak 
stationary gravitational field [21]. We suppose that the metric is close to the locally 
Minkowskian metric, 

g^(x, y) = v»»(y) + v)i 

(14) 

ViwW = f(y) * diag{l, -1, -1, -1} , 

where \h^ u \ <C 1. 

The field is stationary and all time derivatives of g^ vanish. The Cartan tensor is 
symmetric in its three indices. To first order in h^, we have 

r i 1 dhpo 1 df O dhoo n - 
00 2f dx* 4Pdy o!J dx* ' { } 



and 



ffic 00 = 5j» (16) 

All particles here are moving slowly. The velocity y° equal to the speed of light ap- 
proximately and can be regarded as constant. Then, the connection T l 00 reduces as 

L °°~2fdx^ {U) 

From the definition of B^ va p, to first order in h^, we know that all components of B 
vanish. By making use of the non-relativistic perfect fluid approximation, we obtain 
the non vanished component of energy- momentum tensor Too = P9oo- Here p is the 
proper energy density. After the above manipulations, the velocity four-vector reduces 
to the normal velocity v, which is defined in the manifold M and function of distance. 
In the above approximations, the field equation([3]) reduces to 



d ( 1 dh 



oo 



4nGpf(v). 



dxi \2f(y) dx i 
The geodesic equation in Finsler space is given as [21 

d 2 x x i dx^ dx v 



At* ' ^HFW = °' (19) 
Though the form of the geodesic equation is the same with the one in Riemann geom- 
etry, but one should notice that the connection here depends both on coordinates 
and velocities. We will show the physical meaning of the velocity dependence in an 
indirect way by solving the geodesic equation in weak field approximation. In the 
approximation of moving slowly and weak field, the geodesic equation (|T9!) reduces to 



0. 



dH_ 
d^ 

d 2 x % 1 dh 00 f dt\ 



dr 2 2f dx i \dr J 



3Z • (20) 
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The solution of first equation in (|20|) is dt/dr = const.. Dividing the second equation 
in (120]) by (dt/dr) 2 , we obtain 

d 2 x i 1 dh Q _ _<hp_ 

dt 2 ~ 2f dx* dx v [ ' 

where <p is the gravitational potential. Substituting the relation fl2TT) into equation 
(fT8l), we get 

VV = 4nGpf(v). (22) 

The velocity v is a function of coordinate. In spherical coordinate, we can rewrite f(v) 
into p a eo(r). The integration of ff22|) gives 

Vfp(r) = ^ y p(r)p 9eo (r)^ = -£i (23) 

One can see clearly that the effective mass M. is different from the baryonic mass. It 
is enlarged by the factor p ge o( r ) due to Finsler geometrical effect. 

We suppose the energy density of baryons is constant for convenience. Denotes M 
the mass of visible matter. 

We would like limit the metric (fHl) to be the form 



/(0 



5 3ap , i 

6 4yrGp£ ± 



3ap I i 

^G P i f ^24) 

|7rGa pe 3 +(|vrGp) 2 e 4 ) 1/4 • 

Here ao is the deformation parameter of Finsler geometry. It is the measurement 
of deviation from Riemann geometry. The deformation of Finsler space may have 
cosmological significance. One wishes naturally the deformation parameter relates 
with the cosmological constant A. In fact, Milgrom observed[7] that 27rao ~ c(-v/A/3). 
Eqs (122|) and ( 1231 give the geometrical factor of the density of baryons, 

5 _|_ GM_ 

P 9 eo(r) = , " r2a ° ■ (25) 

1 'GM_\2 _i_ GM 
• r 2 ao ' r 2 ao 

In the zero limit of the deformation parameter, familiar results on Riemann geometry 
are recovered. Thus, the acceleration a is expressed as 



GM r 2 + GM/a 
V ^ (r) "V GM j a • (26) 



We can rewrite it in a compact form, 



GM (GM\ 
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where 



v(x) = ./i±£ „ ( 1 *» * » ; . (28 ) 

w V x \ X" 1 / 2 for x < 1 v 7 

It is obvious that the above formula is just the MOND [25]. The MOND reduces to 
Newton's gravity when x 3> 1. Meanwhile, in our model the metric 77^ returns to 
Minkowskian metric. 

Persic, Salucci and Stel analyzed about 1100 rotational curves of spiral galaxies [26]. 
They gave a universal formula for the rotation curves, called Universal Rotation Curves 
(URC), 



VuRc(r) = V{r opt ) 



0.72 + 0.44 log' 



L \ 1.97a; 



1.22 



L* J (x 2 + 0.78 2 ) 1 - 43 



+ L6e -0.4(L/£.) 



x 2 + 1.5 2 (L/L*)°- 4 



1/2 

km/s, (29) 



where r opt is the radius encircling 83% of the integrated light and x = r/r opt . The lu- 
minous matter in spiral galaxies is composed by spheroidal bugle and an extended 
thin exponential disk. The spheroidal bugle has density distribution of the form 
I(r) = J e -7 ' 67 ( r / re ) 1/4 (where r e is the half-light radius) [27] . The thin exponential disk 
has density distribution of the form I(r) = Ioe~ ra (a~ 1 is the disk scale- length) [28] . 
Phenomenological similarities between the URC and MOND have been investigated in 
detail by looking for properties predicted in one framework that are also reproducible 
in the other one [29]. We wish a general Finslerian metric be found to reproduced the 
URC naturally. In fact, many works about dark matter have been done, and the den- 
sity profile of dark matter halos surrounding galaxies have been given. The Finsler 
geometrical factor is determined by 

p geo = 1 + Pdm/p- (30) 
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